Extracting and Stabilizing the Unstable State of Hysteresis Loop 
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A novel perturbation method for the stabilization of unstable intermediate states of hysteresis 
loop (i.e. S-shaped curve) is proposed. This method only needs output signals of the system to 
construct the perturbation form without delay-coordinate embedding technique, it is more practical 
for real-world systems. Stabilizing and tracking the unstable intermediate branch are demonstrated 
through the examples of a bistable laser system and delay feedback system. All the numerical results 
are obtained by simulating each of the real experimential conditions. 
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I. INTRODUCTION 

Bistability [Q-Q has attracted much attention not only 
in the experimental systems but also in the theoretical 
research in the last twenty years for its possible applica- 
tion as an optical device and switcher and so on. As a 
model problem in non-equilibrium statistical mechanics 
bistability has also been considered as an interest- 
ing subject. Usually, we get an S-shaped solution curve 
(i.e. hysteresis loop) in bistable systems, where the up- 
per and lower branches are stable and the intermediate 
branch is unstable. When the mathematical model of 
the system is known, the values of the intermediate may 
be calculated by certain approaches. However, for the 
experimental systems in the absence of a priori mathe- 
matical system model, getting the accurate positions of 
the unstable intermediate state is experimentally inac- 
cessible under ordinary circumstances. How to get the 
unstable intermediate state is considered as a particu- 
larly interesting subject recently |^-^. 

Zimmermann et al. Q firstly obtained an unstable in- 
termediate branch of hysteresis loop (HL) in a thermo- 
chemical experiment about the gaseous reaction mixture 
S2OQF2 ^ 2SO3F, even though they only achieved a 
part of the intermediate branch in the experiment. Their 
method is to introduce an external delay feedback for a 
system parameter to stabilize the unstable intermediate 
branch. As they demonstrated, their method appeared 
suitable for only short delays; but for longer delays in the 
external feedback, some periodic attractors (not present 
in the original system) were predicted to exist by this 
method. So in fact this external delay feedback control 
technique changes the original system dramatically. 

In 1990, Ott, Grebogi and Yorke (OGY) [§ proposed 
a new method of controlling a chaotic dynamical sys- 



tem by stabilizing one of the many unstable orbit em- 
bedded in a chaotic attractor, only through small time- 
dependent perturbations in some accessible system pa- 
rameters. This method was extended by Romeiras et al. 
pO| . For an n-dimensional system 



they used a linear feedback 

Pi+i =po + k(z,;+i - z*) 



(1) 



(2) 



to stabilize the target orbit z,, where k, z e R", and 
p is one of the adjustable system parameters. OGY 
method has been proved unable to change the original 
system. Later researchers extended this method to con- 
trolling the unstable state by using time delay coordi- 
nates technique Principally, OGY method can also 
stabilize the unstable intermediate branch of hysteresis 
loop (HL), but a potential difficulty with a nonadaptive 
application of this technique to real- world system is that 
the unstable state must be known before the use of this 
method. For experimental, especially "black-boxes" sys- 
tems (i.e. little or no information exists regarding their 
governing equation), it is very difficult to determine the 
unstable state exactly. When the OGY feedback (0) 
invalidates in the experiment systems, Z. Gills et al. |l2] 
used p{t) = po + k[z(t) — Zref], ( Zref IS a "guess" value 
for the control objective) and S. Bielawski et al. used 
p{t) = Pa + fiz in Ref. [n3| respectively, to stabihze the 
unstable states. Some of us addressed another method 
based on the OGY method in the variable-parameter 
space to control chaos in Ref. iQ. We suggested a feed- 
back form 



Pi+l =P0+ k(z,;+i - Zi) + k'{pi - Po), 



(3) 
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where k and k' are control coefficients. It is different from 
the previous tracking techniques [p^ , p^ since it need not 
"guess" or "predict" the control target and the unsta- 
ble state in the system can be stabilized automatically 
by choosing suitable control coefficients. 

In 1993, Macke et al. Q used a continuous feedback 
method to acquire the unstable intermediate state of HL 
in an optical bistable system. As a challenge for the 
high-dimensional systems (spatiotemporal systems), Hu 
and He |^ numerically stabilized unstable state of bista- 
bility in a spatiotemporal system described by a par- 
tial differential equation with infecting negative feedback. 
Recently, we successfully expanded our previous method 
[Q to stabilize the unstable state existing in a bistable 
system p6[ | and simultaneously pointed out that there 
exists a transformational law about control coefficients 
between the stable and unstable states by using feedback 
form (^. Based on the Ref. here we show a more 
applicable method to control the unstable intermediate 
state in bistable systems. 

This paper is organized as follows: In Sec. II, the ba- 
sic idea of the control method is shown and some advan- 
tages of the technique are also presented. In Sec. Ill, we 
demonstrate the effectiveness of our technique by apply- 
ing it to five- and infinite-dimensional bistable systems 
respectively. At the end of the paper, a summary of this 
paper and the prospect of applications for experimental 
systems are discussed. 



II. THE BASIC IDEA OF THE CONTROL 
METHOD AND ITS ADVANTAGES 

Most experimental systems are essentially "black- 
boxes" for which one can only get a series of output sig- 
nals zi, Z2, Z3, and measure some adjustable system 
parameters, here z^ £ R". Some researchers have sta- 
bilized unstable states in some experimental systems by 
using a delay-coordinate embedding technique . Here 
we illustrate a well-suited method for experimental sys- 
tems to control the unstable steady states without this 
embedding technique. 

To simplify the analysis we first consider a one- 
dimensional time series, which can be described by a dis- 
crete dynamical system as follows 

Xi+i = f{xi,p,c), (4) 

where x £ R, and / is unknown but is sufficiently smooth. 
Here p is considered as the parameter that is available for 
external adjustment but is restricted to lie in some small 
interval, — Poll < (5, around a nominal value po- c is 
the bifurcation parameter (maybe p and c are the same 
parameter). Let a;*(po) be a fixed point of (i.e., con- 
trolling objective). For the one-dimensional system, the 
feedback form (^) can be chosen as 

Pi+i = Po + k{xi+i - Xi) + k' {pi ~ po) (5) 



to stabilize and track x,, where k and k' are control co- 
efficients. 

Hysteresis loop is a typical bistable phenomenon as 
sketched in Fig 1(a). In this figure, Xs denotes the stable 
lower and upper branches and Xu is the unstable inter- 
mediate branch. The two turning points of the stable 
and unstable states are tangential bifurcation points. Cb 
denotes the value of bifurcation parameter at the tangen- 
tial bifurcation point and Xb represents the corresponding 
value of a variable, ci is a neighboring parameter value 
of tangential bifurcation point, denotes the value of 
variable located in the stable lower branch at c — ci and 
is the corresponding value located unstable interme- 
diate branch. For a hysteresis loop, we have developed a 
method in Ref. to stabilize the unstable intermediate 
state by using the characteristic of the tangent bifurca- 
tion. When applying the feedback law (||) to a bistable 
system, there exists a transformational law between the 
stability region of Xci and the stability region of x'^_^ in 
the control coefficients space (for detail see Ref. Q). We 
have proved that the two regions have no intersection, 
as sketched in Fig. 1(c). One can target the unstable 
intermediate branch of HL according to the method in 
Ref. from which we derived a transformational law 
turning the coefficients obtained at the stable state into 
those suitable for stabiHzing the unstable state, but it is 
complex because a matrix must be calculated to get the 
control coefficients k and k' of stable state. To avoid the 
difficulty, here we present a very easy method to get the 
suitable k and k' without calculation. 

Firstly, to determine k\ we have proved the fact in 
Ref. jl6| that the two stability regions of a;^ and x'^_^ are 
separated by a supercurve k' = 1 and the stability region 
of Xci is located below the supercurve while that of x'^_^ 
is located above when c = ci [see Fig. 1(c)]. k' = 1 — S 
is located in the stability region of Xci , and k' — 1 + 6 
lies in the region of x'^_^ , where 5 is a small positive num- 
ber. If one applies fc' = 1 — i5 to the system, is stable 
while x'^_^ is unstable. On the contrary when one uses 
k' = 1 + 5 to the system, becomes unstable and x'^_^ 
becomes stable. So we can easily choose k' = 1 ± S to 
stabilize the unstable or stable state. 

Secondly, in order to determine fc, one should accept 
the fact that the point (fc, fc') = (0, 0) in the control coef- 
ficient space must be located within the stability region 
of Xci, as sketched in Figs. 1(c) and (d). It is easy to 
understand that because in this case Pi+i is po according 
to the feedback form (||), the system must stabilize the 
Xci state naturally. At the tangential bifurcation point 
Cf), the origin of coordinate is located at the boundary of 
the region [Fig. 1(b)]. Thus one can reasonably expect 
that the stability region of must lie in the neighbor- 
hood of the point (fc, fc') = (0, 0) and can be easily found 
by scanning the coefficient space around the origin point 
(0,0), [see the down triangle in Fig. 1(c)]. One can take 
the values of fc and fc' in the region of , say A(fci, 1 — 15), 
where fci should be chosen in the middle value of region 



2 



in the direction of fc at fc' = 1 — (5 in order to assure that 
A'{ki, 1 + S) is also located in the stability region of x^^. 
If (fc, k') — (fci, 1 + S) is used to replace (fci, 1 — S), the 
perturbations can drive the system dynamical tra- 
jectories to depart from the stable lower state Xc-i and 
terminate on the unstable intermediate one x'^_^ at ci au- 
tomatically. In this case the lower branch loses stability 
and the unstable intermediate state becomes stable. Tak- 
ing {k,k') = {ki,l + 5), one can track the intermediate 
branch when c decreases. 

Finally, if the control coefficients {k,k') — (fci, 1 + S) 
lose control ability when bifurcation parameter c de- 
creases at c = C2 [see Fig. 1(d)], one can use ^2 to replace 
the control coefficient fci, and apply it to the system ac- 
cording to the feedback law (||) to stabilize the unstable 
intermediate state continuously. The reason is demon- 
strated in the following. As sketched in Figs. 1(c) and 
(d) , with the bifurcation parameter c being far away from 
the tangential bifurcation point Cb, the stability regions 
of Xs and also moving along in the positive or 

negative direction of k (for different systems, the mov- 
ing direction is different. Without loss of generality, we 
suppose the region moves to the negative direction of fc 
in Fig. 1). At the same time, the size of regions is also 
changed, but the origin of coordinates (0, 0) is always lo- 
cated in the stability region of Xs (the reason has been 
demonstrated in the last paragraph). At c = C2, though 
A{ki, 1 — (5) is still located in the stability region of x^ 
(the down triangle), A'{ki, 1+6) does not lie in the stabil- 
ity region of Xc2 (the up triangle), so A{ki, 1 + S) is out of 
control for the intermediate state. However, B'{k2, 1 + 5) 
from the corresponding point B{k2, 1 — S) is located in the 
stability region of Xc2 , thus it can continuously stabilize 
and track the intermediate branch. This step can be re- 
peated with the further decrease of parameter c until the 
whole intermediate branch of the HL is obtained. Gen- 
erally (in the following applications of this paper), only 
by choosing fei and fc' at ci , one can get the whole inter- 
mediate state of hysteresis loop by this technique. This 
just verifies the robustness of our method (for detail see 
Sec. IIIA). 

For an n-dimensional system z,;+i = f{zi,p,c) with 
several unstable directions, where z e R" , one can use 

Pi+i =Pq+ k(zi+i - Zi) + k'{pi ~ po) (6) 

to control the unstable intermediate state of HL, here 
k £ R". We can choose ki{ki G k, i = 1,2,3, ...ti) ac- 
cording to the method of choosing fc in a one-dimensional 
case, fc' can still be chosen as 1 — (5 for the stable state 
and 1+6 for the unstable one. In the examples of Sec. 
HI, we consider two dynamical systems which instabil- 
ity is characterized by only one unstable direction. Thus 
the stabilization of the unstable intermediate state of HL 
only requires one system variable to be monitored and 
one system parameter perturbed, so only two control co- 
efficients fc and fc' are used, as done in a one dimensional 
case, to realize the control. 



Usually, for the systems whose mathematical model 
are known, the unstable intermediate steady state of HL 
can also be achieved numerically by Newton's method, 
with the help of the priori knowledge about state func- 
tions suitable to the dynamical system. But for the ex- 
perimental systems whose governing equations unknown, 
there is usually no way to know the state functions, so 
the Newton's method will invalidate. Nevertheless, the 
approach in the present paper whose goal is to provide 
a technique for the experimental systems to stabilize the 
unstable intermediate state of hysteresis loop, stabilizes 
the unstable state by making small perturbations to an 
accessible system parameter, thus it is well-suited for the 
real-world experimental system. 

In addition, this method has other advantages: (i) 
It can leave the original system unchanged because it 
is an extended version of OGY method and inherits its 
characteristic, (ii) The feedback form (^) develops OGY 
method since one does not need to know ( "guess" or "pre- 
dict") the desired orbit x* (which, in fact, is difficult to 
obtain in experiments), and what one needs only to do is 
to record output signals of two times in succession (i.e. 
Xi, Xi+i) and the parameter pi in this step. Thus one can 
obtain the parameter value Pi+i for the next time to con- 
trol the system according to the feedback form It is 
not difficult to design the corresponding experiment ele- 
ments based on the original system, where fc and fc' can be 
considered as adjustable signal amplifiers. When the con- 
trol coefficients fc and fc' are chosen suitably, the system 
will be stabilized onto the control target by itself, (iii) 
Using the outputs of the system to construct the control 
law without the delay-coordinate embedding technique, 
one does not need to know the characteristic of tangent 
space and does not need to calculate the derivative ma- 
trix of Eq. (^) to get the control coefficients. This gives 
us a convenient way in the experimental system in which 
mathematical model is absent, (iv) This technique has 
good resistance to noise. Some experimental systems are 
influenced by noise (e.g. some system parameters change 
as a function of time influenced by temperature, etc.), 
thus the position of the desired orbit moves a small 
distance due to the change of system parameters. Since 
Eq. (|^) does not include x, explicitly, the perturbations 
with any point in the stability region can force the sys- 
tem trajectories to evolve towards the exact position of 
X* (interested readers are referred to Refs. ]l|,|l§). 



III. APPLICATION IN BISTABLE SYSTEMS 

A. Five-dimensional model: An optical bistable 
system. 

Now, let us consider an optical cavity flUed with a 
passive medium, containing of homogeneously broadened 
two-level atoms, and driven by an external coherent laser 
held 1,1,0- Considering only the single mode case, 
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applying the plane wave approximation and taking the 
mean- field limit, one can reduce the Maxwell-Bloch equa- 
tion to 

X = -k[{l + ie)x -y + 2Cp, 
p^xD-{l + iA)p, (7) 
D ^ ~j[{x*p + xp*)/2 + D - 1], 

where the variables x,p are the complex output field 
and the atomic polarization. D is the real population 
difference. Thus the equations are essentially five di- 
mensional. The system parameters C, k, and 7 respec- 
tively represent the small-signal atomic gain, the cavity 
linewidth, and the atomic population decay rate. Both 
K and 7 are scaled to the homogeneous linewidth 7_l [i.e. 
k = (k/7_l),7 = (7||/7_l)]. Given the frequencies of the 
external field, the atoms, and the cavity as too,^^, and 
zuc, we scale the atomic detuning and the cavity mis- 
tuning as A = (wa — ■n7o)/7± and 9 — (nj^ — Wq)/k. 
The time t is measured in units of polarization relax- 
ation time, 7j^. y is the amplitude of the external field, 
an adjustable system parameter. 

Under appropriate conditions, the system exhibits a 
bistable steady state behavior, and the transmitted light 
varying discontinuously with the incident field shows a 
hysteresis loop Using the fourth-order Runge-Kutta 
integration and with fixed step size dt = 0.01, we get a se- 
ries of output signals about the output field x from Model 
. We numerically plotted the HL of this system gov- 
erned by Eq. (0) at C = 30, 6* = 20, 7 = 2, A = 4, and 
k = 0.5, as illustrated in Fig. 2(a). These parameters 
can be realized in the experiment of Ref. [Q. Clearly, 
when the bifurcation parameter y increases from zero to 
11, the real part of output variable x (i.e. Hex) jumps 
up to the upper steady state at y = 9.43. Once on the 
upper branch, if the bifurcation parameter is slowly de- 
creased to zero, Rex falls down to the lower steady state 
at y = 7.0. Conventional stability analysis finds the up- 
per and the lower branches of the S-shaped curve to be 
stable. The intermediate branch is found to be unstable 
and is experimentally inaccessible under ordinary circum- 
stances. When the adjustable parameter y — 9.42, which 
lies in the neighborhood of the turning point from the 
lower to the upper state, we can easily get the stability 
region k e (—28.4, —8.9) at k' = 0.98 by using control law 
yi+i = yo + k{Rexi+i-Rexi) + k'{yi - yo) in the system. 
In our work, we use = — 16 which lies in the middle of 
the stability interval, apply k' — 1.02 to replace k' — 0.98 
to stabilize the unstable intermediate state, and we find 
the system stabilizes automatically to the (original un- 
stable) intermediate state of HL. Figure 2(b) shows the 
targeting process versus integral times n. At the begin- 
ning of the process, Rex is located in the stable lower 
state without control; at n = 12000, we begin to add a 
control with k = — 16, fc' — 1.02 to the system, we find 
Rex increases slowly and is completely located in the in- 
termediate (original unstable) state at n = 19000; when 
n = 40000, we use fc' = 0.98 to replace 1.02 and fc is 



still —16, we observe that Rex decreases and turns back 
to the stable lower state at n = 60000. In fact, when 
the system lies in the intermediate state, we also cancel 
the control and wait the system to restore to the lower 
state slowly. When system stabilizes the intermediate 
state at y = 9.42, then decreasing the bifurcation pa- 
rameter y from 9.42 to 7.0 and using the control law 
with k — — 16, fc' = 1.02 continuously, one can stabilize 
the whole unstable intermediate state of HL, as shown 
in Fig. 2(a). Without modifying the control coefficients, 
we can stabilize and track the whole intermediate state. 
Obviously, this technique is robust for the experimental 
systems. 



B. Infinite dimensional model: A delay feedback 
differential system. 

The method above can also be used in the experimen- 
tal system designed by Gibbs and Hopf et al. (see Fig. 
1 of Ref. [Q). The delay- line hybrid optical bistable 
device can be described by a delay differential equation 
with an infinite number of degrees of freedom, 

T'ip{t) = -ip{t) + A^{1 + 2Bcos[ip{t - tfl) - ipo]} (8) 

where ip = ttV/Vh,V is the voltage fed to the modula- 
tor, Vh is the half-wave voltage of the modulator, r' is 
the composite detector-feedback-PLZT response time, B 
is a coefficient that measures the ability of the modula- 
tor to achieve extinction between the crossed polarizer, 
ipo = — 7r/2 can be set in the experiment, the bifurcation 
parameter = 7r/i = CGiG2lin is proportional to the 
gains Gi and G2 of the amplifiers and to the input laser 
intensity /i„, which is the adjustable parameter in their 
experiment, /i can be measured directly by breaking the 
becking the feedback loop between the final amplifier and 
the modulator. Let t = tnt, (p(t) = ip{t) to rescale the 
Eq. (^, we then have 

r' ■ 

— <^ (t) = -ipit) +Tr^i{l + 2Bcos[if{t - 1) - ipo]}- 

tR 

Denoting r = r'/t/j and dropping the bars from ^, we 
get 

Tip{t) = ~ip{t) +TTfi{l + 2B cos[ip{t - 1) - ^0]}, (9) 

where B = 0.3, as be chosen in Ref. jlSt . 

Equation (^ can be solved by using Adams interpola- 
tion numerically. In order to trace the course of evolu- 
tion, one can draw the evolution curve of the variable ip{t) 
via the time t. We plot the hysteresis loop at t = 6.25, 
/i from to 3.0 and then back to 0, as shown in the Fig. 
3(a). This is under the corresponding experimental con- 
dition of Gibbs et al. [i.e. ta = 160/is t' — 1ms 
in Fig. 1(b) of Ref. (Ts)]. In Fig. 3(a), ^p^ is obtained 
by measuring ip{t) in every time unity. In the increase 
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of /z, the lower branch loses stability at /i = 1.42 and 
the system jumps up to the upper branch. In the oppo- 
site direction, the upper branch loses stability at /i = 0.9 
and the system falls down to the lower branch. There 
exists a bistable region from fj, = 0.9 to ^ = 1.42. Ob- 
viously, our numerical result is in excellent qualitative 
agreement with their experiment. However, the interme- 
diate branch of hysteresis loop [dotted line in the Fig. 
3(a)] has not been observed in their experiment. If the 
control law /i^+i = /uo + k{ip,+i - ipi) + k'{tJLi - fj-o) is 
added to the system (we add a control every time delay 
unit) at Ho — 1.40, here k = —0.6, A:' = 1.02, the lower 
stable branch will be found to lose stability and the sys- 
tem will be locked to the intermediate unstable branch 
ultimately. At this time, if the control is cancelled or 
k = —0.6, k' = 0.98 used to renew the control coeffi- 
cients, the intermediate state is observed to lose stability 
and the system turns in the end back to the lower branch 
via the time t, as shown in Fig. 3(b). When the system 
stabilizes the intermediate branch by fc = —0.6, k' — 1.02 
at /io = 1.40, if one decreases the bifurcation parameter 
fi from /io = 1-40 to /xq = 1-0, one will get the whole un- 
stable intermediate state of the HL, as presented in Fig. 
3(a). 

We also achieve the whole unstable intermediate state 
successfully in another parameter, shown in Fig. 4. From 
this figure, one can find that the upper branch of the 
bistable region has been in the self-pulsing region, but 
we can also scan the intermediate unstable branch by 
the control method in this paper. 



IV. CONCLUSION AND DISCUSSION 

In summary, we have numerically achieved the unsta- 
ble intermediate states of hysteresis loop in two bistable 
systems by a practical control method. Using the feed- 
back law p^+i = pQ + k{xi-f.i—Xi) + k'{pi—pQ) in a bistable 
system with hysteresis loop, the stability regions of sta- 
ble state Xs and unstable intermediate state a;„ have no 
intersection ink — k' plane ]l6t . The two regions are sep- 
arated by a supercurve fc' = 1 and the stability region 
of Xs is located below while that of a;„ is located above. 
For different systems, the stability interval of k is dif- 
ferent but which can be easily obtained by scanning the 
neighbor of zero while the value of k' is same and can be 
chosen as 1 -I- J for the unstable states. One of the main 
advantages of this targeting method is that as the target- 
ing procedure proceeds the stable state becomes unsta- 
ble while the unstable one (i.e., the targeting objective) 
becomes stable and therefore one need not "guess" (or 
"predict" ) the position of the unstable state in advance, 
which certainly benefits the experimentalists. 

On the other hand, this method need not reconstruct 
the system by delay-coordinate embedding technique or 
know the characteristic of the tangential space of the sys- 
tem, and what one only need do is to measure the output 



signals. From output signals, one can construct the con- 
trol law easily. Because of this simplicity, the approach is 
of practical interest for the stabilization of intermediate 
branch in bistable experiment systems. For the experi- 
ments, the key point is to design a group of experiment 
elements to complete perturbation for the adjustable pa- 
rameter Pi+i = po + k{xij^i — Xi) + k'{pi — po) for the 
system. We should like to emphasize that this method 
may be applied by experimentalists to realize the unsta- 
ble output related to bistability from the experimental 
device even if the global mathematical model of the sys- 
tem is unknown. The authors of this paper expect that it 
will be verified in optical or other bistable systems such 
as electronic circle, chemical system, magnetic research 
field, and others. 
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FIG. 1. (a). Schematic of hysteresis loop. Xs denotes the stable branch and Xu denotes the unstable intermediate branch, c 
is a bifurcation parameter and Cb is the tangential bifurcation point, (b)-(d) are the schematices of the stability regions of Xs 
(the solid-line down triangle) and the stability regions of Xu (the solid-line up triangle) at different values of c in the k-k' plane. 



FIG. 2. Controlling of the unstable intermediate branch, (a). Hysteresis loop of an optical bistable system and the controlling 
result of the inermidiate branch. Bar scatters show the lower and upper stable states, and the dotted line denotes the 
intormodiato unstable state, (b) Hex versus integral times n for a representative control technique trial. System switches 
between the lower stable and the intermediate states by using = — 16, fe' = 1 ± 0.02. 



FIG. 3. (a) Hysteresis loop at r = 6.25. Bar scatters show the lower and upper stable states, and the dotted line denotes the 
intermediate unstable state, (b) ip{t) versus time for a representative control techenique trial. The system switches diagram 
between the lower stable and intermediate unstable states by using k = —0.6, fe' = 1 ± 0.02. 



FIG. 4. Hysteresis loop at r = 0.2. The bar scatters and dotted line as same as Fig. 3(a). 
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